! ! tt3f ;?■ 



fa 5 !!!' 



fl-J 

h< 



" l" 



13 P^PCTZPTO 2001 

^09/ 8 0 6 3 2 

WO 00/20632 PCT/GB99/03267 




SUSCEPTIBILITY LOCUS FOR OSTEOARTHRITIS 



Technical Field 

This invention relates to the identification of chromosomal 
5 regions linked to susceptibility to osteoarthritis using 
linkage and association analysis. 

Background of Invention 

Osteoarthritis (OA) is a debilitating disease involving 
10 degeneration of the articular cartilage of synovial joints 5 ' 6 
Although OA has long been considered an inevitable 
consequence of ageing , it has become increasingly apparent 
that OA does have a genetic component . Early-onset forms of 
the disease are associated with several 

osteochondrodysplasias , rare diseases involving abnormal bone 
and cartilage development that are transmitted as Mendelian 
traits 7 . However, the OA in these conditions is secondary to 
the main dysplastic phenotype . The common late-onset form of 
I'm the disease (idiopathic OA) often has no obvious 
RiO environmental (i.e. injury) or characteristic physical cause 
and does not demonstrate a clear mode of inheritance. 
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Over the last 10 years, many genes for single gene or 
monogenic diseases, which are relatively rare in the 

25 population, have been positioned by linkage analysis in 

families, and localised to a small enough region to allow 
identification of the gene. The latter sublocalisation and 
fine mapping can be carried out in single gene rare diseases 
because recombinations within families define the boundaries 

30 of the minimal interval beyond any doubt. In contrast, in 
common diseases such as osteoarthritis, diabetes or asthma 
the presence of the disease mutation does not always coincide 
with the development of the disease: disease susceptibility 
mutations in common disorders provide risk of developing of 

35 the disease, and this risk is usually much less than 100%. 
Hence, susceptibility genes in common diseases cannot be 
localised using recombination events within families, unless 
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tens of thousands of families are available to fine map the 
locus. Because collections of this size are impractical, 
investigators are contemplating the use of association 
mapping, which relies on historical recombination events 
during the history of the population from which the families 
came from. 



Association mapping has been used in over a dozen examples of 
rare single gene traits, and particularly in genetically 
isolated populations such as Finland to fine map disease 
mutations. Nevertheless, association mapping is 
fundamentally different from straightforward linkage mapping 
because even though the degree of association between two 
markers or a marker and a disease mutation is proportional to 
the physical distance along the chromosome this relationship 
can be unpredictable because it is dependent on the allele 
frequencies of the markers, the history of the population and 
the age and number of mutations at the disease locus. For 
rare, highly penetrant single gene diseases there is usually 
one major founder chromosome in the population under study, 
making it relatively feasible to locate an interval that is 
smaller than one that can be defined by standard 
recombination events within living families . The resolution 
of this method in monogenic diseases in which there is one 
main founder chromosome is certainly less than 2cM, and in 
certain examples the resolution is down to 10 0 kb of DHA 
(Hastbacka et al . (1994) Cell 78,1-20). 



In common diseases like OA, which are caused by a number of 
genes or polygenes acting together in concert the population 
frequency of the disease allele may be very high, perhaps 
exceeding 50%, and there are likely to be several founder 
chromosomes, all of which impart risk, and not a 100% 
certainty of disease development. Because association 
mapping is dependent on unpredictable parameters, and because 
founder chromosomes will be several and common in frequency 
in the general population, the task of fine mapping polygenes 
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is currently one of some controversy, and many doubt the 
feasibility at all of a systematic genetic approach using a 
combination of linkage and association mapping. Recently, 
Risch and Marakandis have provided some mathematical 
5 background to the feasibility of association mapping in 

complex diseases (Science 273 1516-1517, 1996) but they did 

not take into account the effect of multiple founder 
chromosomes . 

10 It has often been noted in epidemiological studies that there 

is a female preponderance for OA 5,6 . This may be accounted 
3 for by differential effects on the two sexes of environmental 
*:f factors. However, a Finnish twin study has suggested that 



analysis 3 . Not only have differences in risk between females 
and males been reported but it has also been suggested that 
there are differences in heritability between joint 
groups 2 ' 10 ' 11 . These differences could be the result of genetic 



Summary of the Invention 

The present inventors have identified a region on chromosome 
llq that may harbour a susceptibility locus for OA. This 

25 region was identified following a two-stage non-parametric 
linkage analysis. The first stage involved a random screen 
of the genome using 272 microsatellite markers in 297 OA 
families. The second stage was more selective and involved 
genotyping an additional 184 families for those markers that 

30 demonstrated moderate evidence of linkage in stage 1. This 

revealed three micro satellites on chromosome 11 for whom the 
evidence for linkage increased as the number of families 
studied increased. Finer mapping in and around these 
microsatellites was performed which provided enhanced 

35 evidence for linkage and enabled linked regions to be 
defined. Stratification analysis suggested that the 
chromosome 11 locus may have differential penetrance between 




genetic susceptibility may be greater in women than men 9 . 
This result has recently been supported by a segregation 




locus heterogeneity between the different joints. 
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males and females and between the two different joint groups 
examined (hips and knees) . 

The chromosome 11 linkage appears to be accounted for 
5 principally by females: the highest single point LOD score 
obtained was 3.32 (p=0. 00005) for marker D11S901 in affected 
female-only pairs. In male -only pairs the LOD score for 
D11S901 was 0.04 (p=0.33). The maximum multipoint LOD (MMLS) 
score in our un-stratif ied data was 3.15 in-between markers 
10 D11S1358 and D11S35 whilst the highest MMLS in our stratified 
data was in our female/hip pairs with a MMLS of 3.03 in- 
between markers D11S1314 and D11S916. Linkage disequilibrium 
analysis identified excess transmission of the 258bp allele 
5 (allele 11 as listed in the GDB) of marker D11S937 (5cM 
p5 proximal to D11S901) in our un-stratif ied data (p<0.00!), in 
m our female-only pairs (0 . 01>p>0 . 001) , in our hip-only pairs 
ill (p<0.001) and in our female/hip pairs (p<0.001) . This data 
indicates that the OA susceptibility locus on llg is located 

•liiiStj. 

j|l close to D11S937 and that the 258bp allele of D11S937 may be 
Mo present on a predisposing ancestral chromosome. 

:t , 

j«& Some evidence has also been obtained for linkage with the 

polymorphic chromosome 6 marker D6S2 73 and the polymorphic X 
chromosome marker DXS1068. This may indicate the presence of 

25 loci on these chromosomes that also have a role in OA 
susceptibility. These loci may interact with the locus 
identified on chromosome 11 and play a role in the 
differential penetration observed between males and females. 

30 Open reading frames (ORFs) are stretches of genetic sequence 
which are candidates for being expressed genes . They can be 
identified by the presence of sequence elements which are 
characteristic of coding sequence, such as sequence elements 
from exon-intron boundaries, transcriptional initiation and 

35 termination motifs and start and stop codons . Because large 
amounts of sequence can be screened rapidly for these 
elements, the identification of ORFs is commonly an initial 
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Microsatellite marker loci are designated in this application 
according to the nomenclature conventional in the field of 
5 Human Genetics. This provides a unique designation which 

specifically and unambiguously identifies each marker locus. 
Mapping data for any particular marker locus is readily 
available from conventional sources, such as the Gopher 
server at the Human Genome Mapping Project Resource Centre 

10 (Host = gopher . homo . mrc . ac . uk ? Port = 70 or URL: 

gopher: //gopher . hgmp .mrc .ac.us: 70/ or by anonymous ftp from 

O f tp . hgmp . mrc ,uk : /Oxford Primers ) - 

The present invention relates to chromosomal regions linked 
fife to genetic sequences which affect susceptibility to 
Si osteoarthritis. 

si 

g ■{*';? 

!L. A first aspect of the present invention is a method for 

Ifl identifying individuals susceptible to osteoarthritis 

*¥o comprising obtaining a sample of genomic DNA and detecting 

n the presence or absence of the 258bp allele of D11S937 from 

H ; chromosome 11 . 

Another aspect of the present invention is a method for 
25 identifying female individuals susceptible to osteoarthritis 
comprising obtaining a sample of genomic DNA and detecting 
the presence or absence of the 258bp allele of D11S937 from 
chromosome 11 . 

30 Another aspect of the present invention is a method for 

identifying individuals susceptible to osteoarthritis of the 
hip comprising obtaining a sample of genomic DNA and 
detecting the presence or absence of the 258bp allele of 
D11S937 from chromosome 11. 

35 

A further aspect of the present invention is a method for 
identifying females susceptible to osteoarthritis of the hip 
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comprising obtaining a sample of genomic DNA and detecting 
the presence or absence of the 258bp allele of D11S937 from 
chromosome 11. 

5 Another aspect of the present invention is a method for 

isolating genetic loci associated with susceptibility to OA 
comprising screening a genomic library from an individual who 
is homozygote for the 258 bp allele of D11S937 and 
identifying open reading frames in regions adjacent to said 
10 allele. 

:! ;;^ Another aspect of the present invention is a method for 

\Q isolating genetic loci associated with susceptibility to OA 

•Jf comprising identifying open reading frames in regions 

iffe adjacent to D11S93 7 and comparing said open reading frames in 

"■&f individuals carrying a 258 bp allele of D11S93 7 with said 

ijp open reading frames in individuals with a non-258 bp allele 

B1 , Of D11S93 7. • 

I If! 

S§ Another aspect of the present invention is a method for 

|K isolating genetic loci associated with susceptibility to OA 

comprising screening genomic libraries with the 25 8 bp allele 
of D11S93 7 and identifying open reading frames located within 
500 Kb of D11S937, or more preferably within 100 Kb of 
25 D11S937, or even more preferably within 50 Kb of D11S93 7 or 
most preferably within 10 Kb of D11S937. 

Another aspect of the present invention is the use of the 
258bp allele of D11S937 as a marker for the identification of 
30 loci influencing susceptibility to OA. 

Another aspect of the present invention is the use of the 258 
bp allele of D11S937 as a marker for the ancestral DNA 
variant conferring enhanced OA susceptibility . 

Another aspect of the present invention is the use of the 258 
bp allele of D11S937 as a marker for mapping associated loci 
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to identify genes associated with OA susceptibility. 



Another aspect of the present invention is the use of the 258 
bp allele of D11S93 7 in investigating loci capable of 
influencing susceptibility to OA . 

Another aspect of the present invention is a method for 
mapping loci which affect susceptibility to OA by comparing 
genomic regions containing the 258 bp allele of D11S937 with 
genomic regions containing other alleles of D11S93 7. 

Another aspect of the present invention is a method for 
determining individual susceptibility to osteoarthritis 
comprising obtaining sample genomic DNA from siblings, at 
least two of which have clinical symptoms of osteoarthritis, 
analysing a region of their genomic DNA comprising the 
polymorphic marker DS11901, identifying allele sharing 
between the siblings as defined by a maximum log of the odds 
(LOD) score of greater than 1 and a p-value of less then 
0.25, and determining individual susceptibility to 
osteoarthritis by reference to the allele sharing. 

Another aspect of the present invention is a method for 
determining individual susceptibility to osteoarthritis 
comprising obtaining sample genomic DNA from siblings, at 
least two of which have clinical symptoms of osteoarthritis, 
analysing a region of their genomic DNA comprising the 
polymorphic marker DS11901, and additionally analysing one or 
more of the following; a genomic region comprising the 
polymorphic marker D6S2 73 and a genomic region comprising the 
polymorphic marker DXS1068, identifying allele sharing 
between the siblings as defined by a maximum log of the odds 
(LOD) score of greater than 1 and a p-value of less then 
0.25, and determining individual susceptibility to 
osteoarthritis by reference to the allele sharing. 



Another aspect of the present invention is a method for 
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determining individual susceptibility to osteoarthritis 
comprising obtaining sample genomic DNA from siblings, at 
least two of which have clinical symptoms of osteoarthritis, 
analysing a region of their genomic DNA comprising a 
5 polymorphic marker and located within 20cM of DS11901, or 

preferably within 15 cM of DS11901, or more preferably within 
10cM of DS11901, even more preferably within 5cM of DS11901 
and most preferably within 2cM of DS11901, identifying allele 
sharing between the siblings as defined by a maximum log of 

10 the odds (LOD) score of greater than 1 and a p- value of less 
then 0.25, and determining individual susceptibility to 

w osteoarthritis by reference to the allele sharing. 



Jif Another aspect of the present invention is a method for 

IfS determining individual susceptibility to osteoarthritis 

?f! comprising obtaining sample genomic DNA from siblings, at 

jjff least two of which have clinical symptoms of osteoarthritis. 

S L analysing a region of their genomic DNA comprising a 

|ff polymorphic marker and located within 20cM of DS11901, or 

if) preferably within 15 cM of DS11901, or more preferably within 

lOcM of DS11901, even more preferably within 5cM of DS11901 
m and most preferably within 2cM of DS11901, and additionally 
analysing one or more of the following; a genomic region 
comprising the polymorphic marker D6S273 and a genomic region 
25 comprising the polymorphic marker DXS1068, identifying allele 
sharing between the siblings as defined by a maximum log of 
the odds (LOD) score of greater than 1 and a p-value of less 
then 0.25, and determining individual susceptibility to 
osteoarthritis by reference to the allele sharing 

3 0 

Another aspect of the present invention is a method for 
determining individual susceptibility to osteoarthritis 
comprising obtaining sample genomic DNA from siblings, at 
least two of which have clinical symptoms of osteoarthritis . 
35 analysing a region of their genomic DNA comprising one or 
more of the polymorphic markers D11S901, D11S903, D11S907, 
identifying allele sharing between the siblings as defined by 
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a maximum log of the odds (LOD) score of greater than l and a 
p-value of less then 0.25, and determining individual 
susceptibility to osteoarthritis by reference to the allele 
sharing . 

Another aspect of the present invention is a method for 
determining individual susceptibility to osteoarthritis 
comprising obtaining sample genomic DNA from siblings, at 
least two of which have clinical symptoms of osteoarthritis, 
analysing a region of their genomic DNA comprising one or 
more of the polymorphic markers D11S901, D11S903, D11S907, 
and additionally analysing one or more of the following; a 
genomic region comprising the polymorphic marker D6S2 73 and a 
genomic region comprising the polymorphic marker DXS1068, 
identifying allele sharing between the siblings as defined by 
a' maximum log of the odds (LOD) score of greater than 1 and a 
p-value of less then 0.25, and determining individual 
susceptibility to osteoarthritis by reference to the allele 

sharing . 

Another aspect of the present invention is a method for 
determining individual susceptibility to osteoarthritis 
comprising obtaining sample genomic DNA from siblings, at 
least two of which have clinical symptoms of osteoarthritis, 
analysing a region of their genomic DNA comprising any one of 
the polymorphic markers; D2S202, D3S1266, D4S231, D4S415, 
D6S260, D6S273, D6S286, D6S281, D7S669, D7S530, D11S907, 
D11S903, D11S901, D17S807, D17S789, DXS1068, identifying 
allele sharing between the siblings as defined by a maximum 
log of the odds (LOD) score of greater than 1 and a p-value 
of less then 0.25, and determining individual susceptibility 
to osteoarthritis by reference to the allele sharing. 
Another aspect of the present invention is a method for 
determining individual susceptibility to osteoarthritis 
comprising obtaining sample genomic DNA from siblings, at 
least two of which have clinical symptoms of osteoarthritis, 
analysing a region of their genomic DNA comprising any one of 
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the polymorphic markers D6S273, DXS1068, identifying allele 
sharing between the siblings as defined by a maximum log of 
the odds (LOD) score of greater than 1 and a p-value of less 
then 0.25, and determining individual susceptibility to 
5 osteoarthritis by reference to the allele sharing. 

Another aspect of the present invention is a method for 
identifying loci conferring susceptibility to osteoarthritis 
comprising screening a genomic library with genetic sequence 

10 derived from one or more of the following polymorphic 

markers; D2S202, D3S1266, D4S231, D4S415, D6S260 , D6S273, 

^ D6S2B6, D6S281, D7S669, D7S530, D11S907, D11S903, D11S901, 
D17S807, D17S789, DXS1068 and identifying those isolated 

j M clones which additionally contain open reading frames 

iB 

US Another aspect of the present invention is a method for 

1% identifying loci conferring susceptibility to osteoarthritis 

4i;; ; ;;■ 

comprising screening a genomic library with genetic sequence 
derived from one or more of the following polymorphic 

Jr ,-rf-y 
t W *1 

W markers; D2S202, D3S1266, D4S231, D4S415, D6S260 , D6S273, 
::f;; D6S286, D6S281, D7S669, D7S530, D11S907, D11S903, D11S901, 

FT +i 

D17S80 7, D17S78 9, DXS1068 and identifying open reading frames 
located within 500 Kb, more preferably within 100 Kb, even 
more preferably within 5 0 Kb or most preferably within 10 Kb 
25 of any of these polymorphic markers. 

Another aspect of the present invention is a method for 
isolating genetic loci associated with susceptibility to OA 
comprising screening a genomic library with sequence derived 
30 from the region between polymorphic markers D11S1358 and 
D11S35 and identifying those isolated clones which 
additionally contain open reading frames. 

Another aspect of the present invention is a method for 
35 identifying loci conferring susceptibility to osteoarthritis 
comprising screening a genomic library with genetic sequence 
derived from one or more of the polymorphic markers D6S273, 
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DXS1068 and identifying open reading frames located within 
500 Kb, more preferably within 100 Kb, even more preferably 
within 50 Kb or most preferably within 10 Kb of either of 
these polymorphic markers. 

Description of Drawincrs 

Figure 1 shows a multipoint log of the odds (LOD) analysis on 
the unstratif ied data . 

Figure 2 shows a multipoint log of the odds (LOD) analysis on 
stratified data from female sibling pairs only. 
Figure 3 shows a multipoint log of the odds (LOD) analysis on 
stratified data from hip pairs only. 

Figure 4 shows a multipoint log of the odds (LOD) analysis on 
stratified data from female/hip pairs only. 

Detailed Description of Invention 
Linkage Analysis 

The present inventors collected 481 families in which at 
least two siblings have undergone joint replacement surgery 
of the hip or knee for severe, end- stage idiopathic OA (Table 
1) . Due to the late onset of the disease, none of these 
families contained parents who could participate in the 
study. In stage l 272 microsatellite markers were genotyped 
in 297 of the 481 families. The microsatellites were 
essentially those used by Reed et al 8 with the replacement of 

certain markers that amplified unreliably. Sixteen markers 
from stage 1 had evidence of linkage at p^O.05 (Table 2). 
These markers were then genotyped in the remaining 184 
families. None of the 16 markers had a p<0.05 in this second 
stage although three had a p<0.10: D2S202 (p=0.07), D11S903 
(p=0.07) and D11S901 (p=0.10) (Table 2). When the data for 
stages 1 and 2 were combined and compared to stage 1 only, 
the combined p-value decreased for 4 of the 16 markers: 
D2S202 (p=0.009 for combined vs. 0*036 for stage 1), D11S907 
(p=0.019 for combined vs. 0.05 for stage 1), D11S903 (p=0.0G4 
for combined vs. 0.017 for stage 1) and D11S901 (p=.0.0003 for 
combined vs. 0.0004 for stage 1) increasing the number of 
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families had therefore increased the evidence for linkage at 
these 4 markers . The present invention is related to the 
chromosome 11 markers (D11S901, D11S907 and D11S903) thus 
identified . 

5 

The three chromosome 11 markers that had p^O.05 in the 
combined data set encompass approximately 55cM of the 
chromosome from llpl3 -llql3 . 2 . Most of the chromosome 11 
markers from stage 1 were not taken through to stage 2 as 
10 they did not reach our linkage criteria of p<0.05. Two of 

these markers immediately flanked the 55cM region whilst one 
was located within it. These three markers were genotyped in 
the 184 families of stage 2. In addition, we genotyped 11 
new markers located within and around this 55cM region in all 
481 families to give us denser coverage of this area of the 
chromosome. Table 3 lists the results for all 17 chromosome 
11 markers that were genotyped in all 481 families. The 
lowest p-value was still for marker D11S901 (0.0003) . A 
multipoint analysis gave a maximum LOD score (MMLS) of 3.15 
distal to D11S901 and in-between markers D11S1358 and D11S35 
(Figure 1) . 
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The results were stratified into six categories : those 
families that were affected female-only pairs (196 families) , 

25 affected male-only pairs (male and/or female) (54 families) , 
affected female- only pairs who had undergone hip replacement 
but not knee replacement (female/hip pairs) (132 families) 
and affected male-only pairs who had undergone hip 
replacement but not knee replacement (male-hip pairs) (72 

30 families) (Table 1) . We did not stratify for female/knee 

pairs or male/knee pairs as the number of families were too 
low (21 and 8 respectively) and thus any significant results 
may simply be the result of stochastic variation. 



35 



Stratification revealed that the linkage to chromosome 11 was 
predominantly accounted for by affected female-only pairs 
with a single point LOD score of 3.32 (p=0 . 00005) for marker 
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D11S901 (Table 4) . Five markers immediately flanking D11S901 
also supported linkage in female-only pairs at p<0 . 05 with 
linkage extending from D11S1314 (p^O.006) to D11S135S 
(p=0.015), a distance of approximately 20cM. The MMLS for 
5 female-only pairs was 2.81 in-between markers D11S901 and 

D11S1342 (Figure 2) . For affected male-only pairs there was 
no evidence of linkage at D11S901 (p=0 . 33) or with any of the 
markers distal or proximal to D11S901 except for marker 
D11S4089 (p=0.043, 40cM distal to D11S901) . There was much 
10 greater evidence for linkage in hip-only pairs at D11S901 

(p=0.0004) than in knee-only pairs (p=0.093) . However, this 
difference may be due to stochastic variation as there were 
only 54 knee-only pairs compared to 311 hip-only pairs. The 
MMLS in hip-only pairs was 2.58, again in-between markers 
lis D11S901 and D11S1342 (Figure 3) . In female/hips the p-value 

y for D11S901 was 0.001 compared to 0.12 in male/hip pairs. 

Til 

I'm The MMLS in female/hips was 3.03 in-between markers D11S1314 

v and D11S916 (Figure 4) . Overall these results suggest that 

f% the linkage to chromosome 11 in our families is restricted to 

fgO females with OA. 

CfK! 

J*& Markers from other chromosomes 

Since stratification analysis of chromosome 11 revealed 
apparent significant differences between the different 

25 categories examined we also stratified the 12 other markers 

that had a p<0.05 in stage 1 and which were not on chromosome 
11. The majority of these 12 markers will of course 
represent false positives from our stage 1 analysis; the p- 
value for each increased when the combined analysis was 

30 compared to stage 1 (Table 2) . Nevertheless, they cannot all 
be discounted. Table 5 lists the stratification results for 
these 12 markers. 

Two of the 12 markers are of interest when one considers the 
35 results of previous studies: D6S273 and DXS1068. 



D6S2 73 maps close to the HLA complex on chromosome 6p and 
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Pattrick et al 1 " have reported association of nodal OA with a 

specific allele of the HLA-A gene. Furthermore, the C0L11A2 
gene, which encodes the a2 chain of the cartilage collagen 
type XI, is tightly linked to the HLA complex. This gene has 
5 been linked to and found to harbour causal mutations for the 
osteochondrodysplasia Stickler syndrome 15 . This syndrome has 
severe, precocious OA as one of its many phenotypic 
components . 

10 Leppavuori et al 13 have reported linkage of DIP OA to 
cytogenetic band Xpll.3. DXS1068 also maps to Xpll.3. 

h D Association Analysis 

K Having established convincing evidence for linkage to 

Sfe chromosome 11, the chromosome 11 markers were tested for 

Iff. association. We used the Transmit software program 16 to test 

01: for association with linkage disequilibrium by the 

1, transmission disequilibrium test. 

U $ 

2:f) For chromosome 11 we analysed the genotyping data without 
f ; ! stratification and with the following stratification 

criteria: female-only pairs, hip-only pairs and female/hip 
pairs. We tested eight markers: D11S901 plus the seven 
markers which immediately flanked this marker: D11S1314, 
25 D11S916, D11S937, D11S1342, D11S1358, D11S35 and D11S898 . 
These eight markers cover the linkage detected in the un- 
stratified and the stratified data sets. 

In the un-stratified data disequilibrium was positive only 
30 for marker D11S937 with a global x 2 4df = 12 .35 (0 . 02>p>0 . 01) 

(Table 6) . Disequilibrium was also positive for this marker 
in affected female-only pairs with global x^df^lS . 19 
(0 . 01>p>0 . 001) , in hip-only pairs with global x z Adf =l%.01 
(0 . 01 . >p>0 . 001) and in female/hip pairs with global x 2 3df= 
35 12.73 (0 . 01>p>0 ♦ 001) . In each case it was the same allele of 
D11S93 7 that demonstrated significant distortion in 
transmission (258bp allele [allele 11] as listed in the GDB) : 
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observed transmission in all families of 254 compared to 
expected transmission of 233.88 (x 2 idf H H * 38 , p<0.001), 
observed transmission in affected female-only pairs of 102 
compared to expected transmission of 91.8 (xSdf = 10 . 01 , 
0 . 01>p>0 . 001) , observed transmission in hip-only pairs of 159 
compared to expected transmission of 139.69 (x 2 idf s 17.74,, 
p<0.001) and observed transmission in female/hip pairs of 65 
compared to expected transmission of 55.8 (xidf s 12.01, 
p<0.001). Marker D11S937 is only 5cM proximal to D11S901, 
the marker which has the strongest single-point evidence for 
linkage in our un-stratif ied and stratified analysis (Table 
4) . Our linkage and association results therefore indicate 
that the OA susceptibility locus on liq is located close to 
D11S937 and that the 258bp allele of D11S937 may be present 
on an ancestral chromosome that harbours a casual DNA 
variant . 



Of the seven other chromosome 11 markers tested, 
disequilibrium was positive in the stratified data for marker 
D11S901, but only in the hip-only pairs with a x%df 15.44 
(0 . 02>p>0 . 01) . This disequilibrium was accounted for by the 
significant (p^O.05) distorted transmission of two alleles of 
D11S901: the 168bp allele (allele 3) and the 170 bp allele 
(allele 5) as listed in the GDB . The observed transmission 
in hip-only pairs of the 168bp allele was 3 84 compared to 
expected transmission of 371.6 (x 2 idf=4.05, 0 . 05>p>0 . 02) , 
whilst the observed transmission of the 170 bp allele was 190 
compared to the expected transmission of 18 0.2 8 (x 2 idf"4 . 10 , 
0.05>p>0.02) . Since two alleles are accounting for the 
positive disequilibrium, it is likely that this 
disequilibrium represents a false positive. 

Methods 

Affected sibling-pairs 

Families were recruited which contained at least two siblings 
two had undergone one or more total hip (THR) and/or total 
knee replacements (TKR) for idiopathic OA. Clinically these 
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patients are at the severe end of the OA spectrum with 
advanced radiological changes. The idiopathic diagnosis was 
supported by clinical, radiological, operative and 
histological findings: patients who had undergone joint 
5 replacement surgery secondary to other factors, such as 
intra-articular fracture or rheumatoid arthritis, were 
excluded. Families were ascertained at three centres within 
the United Kingdom: The Nuffield Orthopaedic Centre in 

Oxford, the Royal Orthopaedic Hospital in Birmingham and 

10 Musgrave Park Hospital in Belfast. Idiopathic OA is 
typically a late-onset disease and parents of affected 
siblings are rarely available. Of the 481 families used in 
the study none contained a parent who was able to 
participate. We therefore collected siblings who had not 
undergone joint replacement to assist in the determination of 
identity~by~descent (IBD) allele transmittance . The 481 
families were comprised of 1052 affected individuals plus an 
additional 302 unaffected siblings (Table 1) . 59.3% of the 
affected individuals were female, 40,7%- were male. For each 

&b individual, 25ml of venous blood was collected into EDTA 

11 ' 

tubes and DNA was extracted by conventional techniques. 



Markers and Qenotyping 

Our initial screening panel of 2 72 microsatellite markers was 
25 essentially the panel used by Reed et al 8 . The additional 

microsatellite markers used to provide a denser coverage of 
chromosome 11 were obtained from the GDB or from the ABI 
Prism Linkage Mapping Set (Version 2) . The markers were 
amplified using conventional conditions with either the 

30 forward or the reverse primer in a PGR pair f luorescently 
labelled. The amplification products were electrophoresed 
through 6% acrylamide using an Applied Biosystems 377 
Automated DNA Sequencer 8 . Alleles were sized using Applied 
Biosystems Genescan version 2.0.2. and Genotyper version l.l 

35 software . 



Linkage and Linkage disequilibrium analysis 
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Initially error checking procedures were employed for all 
families for each marker. After identification of 
straightforward mis - inheritances , more subtle transmission 
errors were detected using the PedCheck program 17 . The entire 
family data set was tested with Relative which is able to 
produce a probability calculation (based on 50 or more 
unlinked markers) that full sibships are in fact half sibs or 
even unrelated (due to unknown adoption of laboratory error) . 
All 481 families used in the study successfully progressed 
through these checks. In addition markers were checked for 
having excess homozygotes , based on their allele frequencies 
and heterozygosities. Markers shown to be unreliable were 
eliminated from the study. Finally the marker data were 
haplotyped for each chromosome using Simwalk2. This checks 
for areas on the chromosome where excessive recombination 
events may alert us to genotyping errors or mis-assignment of 
a marker position. 

Non parametric linkage analysis was performed using the 
SIBPAIR module of the ANALYZE package 18 . This module is able 
to use data from siblings to determine identity-by-descent 
(IBD) allele transmittance . In this way it is appropriate to 
our study since we were unable to collect parents of our 
affected siblings. In the linkage analysis, siblings who had 
not undergone joint replacement were given a clinical status 
of unknown. The SIBPAIR module produces a singlepoint LOD 
score and its asymptotic P-value. Allele frequencies were 
calculated from the input data using either GAS or Downfreq 
(part of the ANALYZE package) . Subsequent multipoint 
analyses was performed using ASPEX which calculates its own 
allele frequencies from the data set, using a maximum 
likelihood method, and employs marker information across the 
chromosome simultanenously 19 . ASPEX produces maximum LOD 
score (MLS) under an additive model. In addition it produces 
an exclusion map along the entire chromosome based on a fixed 
value for X<? . 
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We tested for linkage disequilibrium by the transmission 
disequilibrium test (TDT) using the TRANSMIT software 
program 16 . Alleles with a frequency <0.1 were pooled. A 
global x 2 statistic was computed for each microsatellite . if 
5 a microsatellite was significant a x 2 statistic was then 
computed for its individual alleles . 



Stratification 

We stratified for sex, for joint replaced (hip or knee) and 
10 for sex combined with joint replaced. 

n For those families with more than two affected siblings and 
k D in which the siblings were not of the same sex, the affected 
SJ sibling of opposite sex to a same pair was given an affected 
HE status of unknown in the linkage analysis. In this way we 

ij ij 

JJt were stratifying for sex whilst not excluding siblings who 

'A "ttF- 

ill could be used to assist in the determination of identity-by- 
L descent (IBD) allele transmittance. 

ftp A hip -only pair were siblings who had each undergone joint 
A replacement of the hip only (mono or bi- lateral) whilst a 
IM ; knee-only pair had undergone joint replacement of the knee 

only (mono or bi-lateral) . If an affected pair was composed 
of one sibling who had undergone joint replacement of one 
25 joint type only (hip or knee) whilst their affected sibling 

had undergone joint replacement of the hip and knee then that 
pair were excluded. For an affected trio, if a pair of the 
siblings had both undergone joint replacement of the same 
joint type only (hip or knee) whilst the third sibling had 
30 undergone both hip and knee replacement, then the concordant 
pair were used in the stratification study whilst the third 
sibling was given an unaffected status in the linkage 
analysis. Again, we were attempting to maximise our 
determination of IBD allele transmittance . 



35 
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Table 1 



a) The families used in stages I and 2 together with the combined total 

figure. Also listed are the pairs concordant for different stratification criteria, b) 

A list of the individuals in the study. 



a) 

Families 



Families 

sibling pair 
sibling trio 
sibling quattro 
other* 



Female only 
pair 
trio 



Stage 1 

297 

265 

23 

7 

2 



132 
121 
II 



Stage 2 

184 

150 

27 

5 

2 



64 
59 
5 



Total 

481 

415 

50 

12 

4 



196 
180 
16 



Male only 


61 


41 


102 


pair 


57 


36 


93 


trio 


3 


5 


8 


other" 


1 


0 


I 


Hip only 


194 


117 


311 


pair 


170 


96 


266 


trio 


18 


17 


35 


quattro 


4 


2 


6 


other* 


2 


2 


4 


Knee only 


34 


20 


54 


pair 


33 


19 


52 


trio 


1 


1 


2 


Female/hip only 


85 


47 


132 


pair 


77 


46 


123 


trio 


8 


1 


9 


Female/knee only 


16 


5 


21 


pair 


16 


5 


21 


Male/hip only 


45 


27 


72 


pair 


41 


25 


66 


trio 


3 


2 


5 


other* 


1 


0 


1 


Male/knee only 


4 


4 


8 


pair 


4 


4 


8 
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b) 

Individuals 





Stage 1 


Stage 2 


Total 


Affected individuals 


641 


41 1 


1052 


Female 


394 


230 


624 


Male 


247 


181 


428 


Hip only 


479 


309 


788 


Knee only 


121 


77 


198 


Hip & Knee 


41 


25 


66 


Female/hip only 


287 


172 


459 


Female/knee only 


77 


42 


1 19 


Female/hip & knee 


30 


16 


46 


Males/hip only 


192 


137 


329 


Males/knee only 


44 


35 


79 


Males/hip & knee 


11 


9 


20 


Additional siblings*" 


211 


91 


302 


Female 


107 


49 


156 


Male 


104 


42 


146 



:> i- 



it/:; ? 



B S i( 



'[ SKJi? 
■t .- 



a Other relative pairs such as cousins, uncles, aunts. 

b Since our families lack parents, siblings who had not undergone joint replacement surgery 
were collected to assist in the determination of parental genotypes. These siblings were given 
an <c unknown" clinical status in the linkage analysis. 
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